function covarianceNum = covariance(a,b);

covarianceNum = (var(a+b) - var(a-b))/4;